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The goal of this paper is to study the effect of entanglement on the running time of a quantum 
computation. Adiabatic quantum computation is suited to this kind of study, since it allows us 
to explicitly calculate the time evolution of the entanglement throughout the calculation. On the 
other hand however, the adiabatic formalism makes it impossible to study the roles of entanglement 
and fidelity separately, which means that results have to be interpreted carefully. We study two 
algorithms: the search algorithm and the Deutsch- Jozsa algorithm. We find some evidence that 
entanglement can be considered a resource in quantum computation. 
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INTRODUCTION 



In the past few years, it has become increasingly clear that entanglement is an important resource in quantum 
information science. It is crucial for quantum teleportation 0], and has been proven to be necessary for an exponential 
speed-up of quantum computation over classical computation 0-0- 01- However, there is some uncertainty over the 
precise role that entanglement plays in general. It seems not to be necessary for improvements less dramatic than 
exponential. Lloyd @ has shown that quadratic speed-up can be achieved with no entanglement whatsoever, using 
single particle states, albeit at the very significant cost of an exponential increase in what he called "overhead", eg. 
the number of particle detectors required. A recent implementation of Grover's search algorithm using only classical 
waves [fj suggests that in this case interference is the relevant factor. Finally, it has been shown by Biham et al |7j 
using specific examples that quantum computers without entanglement can still provide improvements that cannot, 
even in principle, be classically achievable . 

One of the complications in analyzing the role of entanglement, as Lloyd's analysis, for example, indicates, is that 
there are many other factors which affect the speed of a quantum computation. Examples are the Hilbert-space 
dimension (e.g. in and even the measurement process itself In the case of adiabatic quantum computation, 
it was shown [23] that a temporary, large input of energy into the system could also increase the computation speed. 
This is similar to the results of Lloyd Q cited above, with recoverable energy playing the role of "overhead". Finally, we 
Q" 1 mention the connection highlighted by Wei et al [Tq| between the minimal running time of an adiabatic computation 
and the fidelity (distance) between the initial and final state in the adiabatic algorithm. 
• i-h , In light of all these results, it is of great interest to analyse as fully as possible the precise role that entanglement 
^ ■ plays in achieving speed-up of quantum algorithms. In the following we will focus on the dynamical role of entangle- 
ment in adiabatic quantum computing. Adiabatic quantum computation provides an alternative model for quantum 
computation which considers a system evolving under a time-dependent Hamiltonian instead of through a sequence of 
gates. It has recently been shown to be equivalent to standard quantum computation with respect to computational 
power (Tlj . and may be more promising when it comes to robustness against decoherence or other types of errors (e.g. 
|12l|). For our purposes it is an ideal theoretical framework in which to study the role of entanglement and fidelity, 
since one can track the actual evolution of these quantities as a continuous function of time. 

Whereas some work has been done studying the entanglement dynamics for "naturally occurring" many-body 
systems, e.g. a Heisenberg chain 0], two spins in a magnetic field or two qubits in a thermal bath 0], the 
behaviour of entanglement in systems while a quantum algorithm is running has only been investigated for a few 
special cases: Orus and Latorre ^| calculated the time dependence of entanglement for two adiabatic quantum 
algorithms (NP-exact cover problem and Grover's search problem). In the former, which runs in polynomial time, 
they found indications that the entanglement grows linearly with the number n of qubits. In the latter case, which 
achieves only a quadratic speed-up, t hey showed analytically that the maximum entanglement is bounded in the limit 
of large n. Finally, Roland and Cerf [Tj| have studied the fidelity as a function of time in the locally adiabatic search 
algorithm. They showed that for the optimized algorithm, which achieves a yN speed-up, the fidelity is linear in the 
optimized time variable, in keeping with the fact that the standard quantum search algorithm changes the fidelity in 
equal increments. 

In addition to shedding light on the specific role of entanglement in the efficiency of quantum algorithms, a systematic 
investigation of the time evolution of entanglement during adiabatic quantum computation may also provide other 
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benefits. Firstly, to achieve a better understanding of the engineering of quantum information processes, it is necessary 
to understand the dynamics of the entangled states which are used for the quantum computation. On a more 
fundamental level, new insight from quantum information theory may lead to a better understandin g of interacting 
many-body systems, and phenomena like quantum phase transitions and critical behaviour (see, e.g., 

We must, at this stage, emphasize an inherent limitation associated with the study of the time evolution of entan- 
glement during an adiabatic quantum computation. We generally start from a set of initial states that have zero initial 
entanglement but different initial fidelity relative to the desired final state. We expect that the initial fidelity (which 
measures the separation of the initial state from the final state) and the amount of entangement that is produced 
during the evolution, will both affect the running time of the algorithm. Because the Hamiltonian, and therefore 
the algorithm, depends explicitly on the choice of initial state, it is not possible to study the role of each variable 
separately. Consequently, it is difficult to determine to what extent each is influencing the other. 

We study two different algorithms, adiabatic quantum search and the adiabatic Deutsch-Jozsa algorithm We 
also discuss a variation of the Deutsch-Jozsa algorithm proposed in [23] which we will refer to as the constant-time- 
Deutsch-Jozsa algorithm. To study the possible dependence of our results on the definition of entanglement, we also 
introduce a new function to describe the physical property of entanglement and compare the results with those for 
the entropy of entanglement. Similarly, we compare the results for the fidelity with those for a different function 
characterizing the distance between states. We also calculate the time evolution of these quantities as a function of 
an optimized time variable, which corresponds to applying the adiabatic approximation in a maximally efficient way 
0- 

For the search algorithm, using the unoptimized time variable, we find that the region where the magnitude of the 
rate of change of the fidelity (as defined in Q48JI) is greatest corresponds to the region of greatest entropy production. 
Using the optimized time variable, the n = 2 initial states that produce more entanglement have longer running 
times. For n = 3 this correlation only holds when we use the new function introduced in this paper to measure 
entanglement. For Deutsch's algorithm, using the optimized time variable, we find that states that generate entropy 
more quickly have a shorter running time. The constant-time-Deutsch-Jozsa algorithm behaves like the search 
algorithm for n = 2: an increased production of entropy is correlated with a longer running time. An earlier form of 
some of these results was described in (2ll |. 

Our paper is organized as follows: In section <|l]} we review adiabatic quantum computation and the algorithms 
that we will study. In section ijHD we review the entropy of entanglement and the fidelity, and introduce two related, 
but new, definitions of functions to describe the same physical properties. The results are presented in section ijlVJI 
and we present our conclusions in section 10- Appendix IpfJ contains some notational definitions. In Appendix {Bj 
we derive expressions for the reduced density matrices. In Appendices 10, 10 and (|Ej we present some details of the 
calculations related to a new definition of entanglement. 



II. ADIABATIC QUANTUM COMPUTATION 

Standard quantum computation uses quantum circuits, in which a sequence of unitary transformations acts on a 
set of qubits. The time evolution takes place in discrete steps. The building blocks of an adiabatic quantum computer 
are qubits, too, but the time evolution of the system is described by a time-dependent Hamiltonian, and is therefore 
much closer to the way a physicist describes interacting (many-body) systems. 

Calculations for a system with a time dependent Hamiltonian are much more involved than for a time-independent 
one, but fortunately several approximation methods exist for different limiting cases. The adiabatic approximation 
applies to systems in which the Hamiltonian varies sufficiently slowly in time. Roughly speaking, this approximation 
says that if a system is in an instantaneous eigenstate and its Hamiltonian evolves slowly enough in time, the system 
will remain in this instantaneous eigenstate. The eigenstate and its energy eigenvalue evolve continuously, but there 
is no transition to a different state. (We have restricted our considerations here to the most simple case, assuming 
that the spectrum of H{t) is entirely discrete and nondegenerate. In the case of a degenerate spectrum, the adiabatic 
approximation generalizes to the statement that each degenerate eigenspace, instead of the individual eigenvectors, 
evolves independently [HI)- This situation is intuitively plausible, since we know that for a time-independent Hamil- 
tonian, the system would stay in a given eigenstate for all times. The approximation is made quantitative by the 
adiabatic theorem, which we briefly sketch below, following the excellent presentation in [22l |. 



3 



A. Adiabatic Theorem 

Consider a system obeying a time-dependent Schrodinger equation 

Am) 



dt 



W)Mt)), (i) 



where here and in the following we set K=l. The time-dependent Hamiltonian has instantaneous eigenstates \E n (t)) 
which satisfy 

H(t)\E n (t)) = E n (t)\E n (t)). (2) 

Note that both the energies and the state vectors are functions of time. To simplify the notation we will suppress the 
T and denote the ground state and its energy by {\E^) , E-}, the first excited state and its energy by , E + }, 

and higher states and energies by {\E n ) , E n }. Using these states as a basis, we can express a general solution of Q 

as 

|*(t))= 5> n (t)e^«|£ n (i)>, (3) 

n 

where 

a n (t) = - f E n {t')d£ (4) 



J o 

is the dynamical phase. After some algebraic manipulations, we obtain the following expression for the time evolution 
of the coefficients: 

a k = -a k (E k \E k ) - V (W < (o.- t ) | (5) 

The condition for adiabatic evolution is that the coefficients a k (t) evolve independently from each other, i.e. that 
their dynamical equations do not couple. This is ensured if the following condition is satisfied: 



max 

0<t<T 



(E k \H\E n 



E n — E k 



< min \E n -E k \ (6) 

0<t<T 



where T is the total evolution time. Physically, this means that if the time evolution of the Hamiltonian (in units of 
the energy gap between the levels under consideration) is slow compared to this energy gap, for all pairs of energy 
levels, the condition for adiabaticity is fulfilled. 

If we consider specifically the adiabatic evolution of an initial state which is an eigenvector, 

|*(0)) = \E m (0)), (7) 

then a m (0) = 1 and a„(0) =0 Vrt ^ m, and the evolution equation reduces to 

a m — —a m (E m \E m ). (8) 

The adiabatic evolution of the state is 

|^( f ))=e i («*»W+T«W)|£; m (i)) J (9) 
where the geometrical phase (or Berry phase) 7 m (0 is given by 

lm(t)=i[ dt'{E m {t')\E m {t')). (10) 
Jo 

This phase depends on the path in parameter space that is traversed by the time-dependent parameters in H. Eqs. 
and © state the content of the adiabatic theorem: An instantaneous eigenstate |£" m (0)) evolves continuously to 
the corresponding eigenstate \E m (t)} at a later time, without any transition to other energy levels. 

When the time dependence of the Hamiltonian is not small enough to fulfill the condition for adiabaticity, transitions 
to other eigenstates become allowed. Equivalently, the dynamical equations © will couple distinct eigenvectors, and 
the eigenvectors will not evolve completely independently. The time evolution of such a system could be obtained by 
solving the system of coupled differential equations 101 . 
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B. The Hamiltonian for adiabatic quantum computation 

Adiabatic quantum computation takes the system from an initial state \ipo) to a final state |V>i), which are eigenstates 
of the Hamiltonians H and Hi, respectively. The initial state is assumed to be easy to construct, for practical 
purposes, and the final state encodes the solution to the computational problem, which is unknown. The time 
dependent Hamiltonian is constructed so that it will drive the system from the initial to the final state, moving 
slowly enough to ensure adiabatic evolution, so that it will end up in the desired state with large probability. This 
Hamiltonian is usually constructed as a linear combination of the initial and final Hamiltonians, 

H(t)=f(t)H +g(t)Hi. (11) 

The functions f(t) and g(t) are arbitrary functions of time, usually chosen to be monotonic, which are subject to the 
boundary conditions 

/(0) = 1; g(0) = (12) 
/(T) = 0; g(T) = l (13) 

where T represents the total computation time. For example, the simplest choice of these functions is 

f(t) = l-s(t); g(t) = s(t); s{t)=t/T (14) 

where we have introduced the dimensionless time variable s(t). The initial and final Hamiltonians are constructed 
from the initial and final states as 

H =I-\^o)(M (15) 
Hi=I-\^i){il>i\ 

so that the ground state of H is the state \ipo) and the ground state of Hi is the state \ipi). 



C. Running time of the algorithm 



We define a small dimensionless parameter corresponding to © : 



max 

0<t<T 



(E-\H\E + ) 



(16) 



For e <C 1, the adiabatic approximation is satisfied. Rewriting lllfit in terms of the dimensionless variable s — t/T 
and assuming that the bound is saturated we obtain: 



T = max 



|<£_|4£|£u 



ds 



o<s<i e{E + - EJf 



(17) 



We obtain the running time by choosing a particular value for e, which represents the accuracy of the calculation. 
Physically, for times greater than T the system can be assumed to be in the final state |V>i) (see Eq. iflU ), up to 
errors of order e 2 . 

The minimal running time can be obtained by optimizing the adiabatic approximation Effectively, we apply 
the adiabatic condition itlfi)) to each infinitesimal time step dt throughout the evolution, instead of globally to the 
entire time interval T. We replace the function s(t) = t/T with a non-linear function whose slope is large in regions 
where the energy gap is large, and small when the energy gap is small. The function s(t) is obtained by solving the 
differential equation: 



dt = ds 



e(E, 



E.) 



(E-\^\E + ) 

We choose e = 0.01 and obtain the optimized running time as: 

(E+-E_f 



Tn 



(0.01) 



ds 



(E-\^\E + ) 



(18) 



(19) 



5 



D. Adiabatic Quantum Search Algorithm 

The goal of the search algorithm is to find a marked item, the "needle," in a large unsorted database, the "haystack," 
in as few steps as possible. Classically, 0(N) steps are needed for a database of size N, since all the objects have to 
be examined, one after another. For a quantum search, the items are represented by the computational basis states, 
one of them being the marked state \m). The unsorted database is represented by an equally weighted superposition 
of these states: 



N-l 

= ^Ew ( 2 °) 







which is motivated by the idea that if a measurement were made on the initial state, the probability of selecting 
any item would be the same. The state |^o) is taken as the initial state for the calculation, and the final state is 
= \m), the marked state that we are searching for. Without loss of generality, we choose \m) — |0). In addition 
to using the conventional "haystack" initial state, we perform the numerical calculation using several other choices. 
Specifically, we vary the initial state in order to study the correlation between running time, entanglement and fidelity. 
For the adiabatic quantum search, the Hamiltonian is constructed from the initial and final state as in Eos. ltTT|l 



quai 

and l(T5)) (see j2^| or [24|]). Its matrix elements are 



Hn = fO--Jf) (21) 
Hu = /(I -4) + 9 farMl 



N' 

N 



f 

= -— fori ^ j 



The lowest two eigenvalues are 




E (') = -( (/ + </)±\/(/-<7) 2 + |f/<? , (22) 



and the only higher energy eigenvalue E n = f + g is (N — 2)-fold degenerate j23|. To calculate the running time, we 
only need to determine explicitly the lowest eigenvector. Within the adiabatic approximation, this lowest eigenvector 
gives the instantaneous state of the system and is used in the determination of the entanglement and fidelity as 
functions of time. 



E. Deutsch-Jozsa Algorithm 

The Deutsch-Jozsa problem j2^| is to determine whether a Boolean function 

F ■ {0, 1}" — * {0, 1} (23) 

is constant (i.e. all outputs are identical) or balanced (i.e. it has an equal number of zeroes and ones as output). 
Classically this requires ~ 0(2™/2) evaluations of the function. Quantum mechanically one can start with an equal 
superposition of all values of the input variable and determine the solution using a single application of the oracle 
operator Uf that is the quantum analogue of the classical function evaluation. Although this is a somewhat artificial 
example, it provides an excellent illustration of how quantum superpositions allow maximal use of the information 
content of the full Hilbert space. What is less clear is whether entanglement plays any role in this case. 

For the simplest case of n = 1, the function would be constant if F(0) = F(l), and balanced if F(0) ^ F(l). The 
four possible outcomes in this case are 

F(0) = F(l) = (constant) (24) 
F(Q) = F(l) = 1 (constant) 
F(0) = 0; F(l) = 1 (balanced) 
F(0) = 1; = (balanced) 

Classically one has to determine both F(0) and F(l) to infer the nature of the function, since the knowledge of one 
does not shed any light on the value of the other. Using the quantum version, Deutsch's algorithm [26], the nature of 
the function can be determined by making one measurement. 
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An adiabatic realization of the Deutsch-Jozsa algorithm has been introduced in The initial state is the same 
as in the search algorithm, Eq. l(20jl . Unlike the case of the search algorithm, there is no natural physical motivation 
for this initial state in Deutsch's algorithm, but it has been shown in 0] that this choice is optimal in the sense 
that the time evolution is independent of the values of a and (defined below) , and thus independent of the result 
of the algorithm. For n qubits, the running time for this algorithm scales as \/~N = 2™/ 2 . Once again, we perform 
the numerical calculation using several other choices for the initial state, in order to study the correlation between 
running time, and entanglement and fidelity. 

The final state is given by 



|Vx> =o|0) + 







N-l 



*:=i 



where 



1 

N 



a~e{0,l}™ 

= 1-a. 



(25) 



(26) 



If F is constant, a — 1 and (3 = 0, and if F is balanced, a = and (3 = 1. This means that if the measurement of the 
final state after the required running time yields |0), the function is constant, otherwise it is balanced. 
The Hamiltonian is given by: 



Ha 



f(N-l) 
N 

f(N-l) 



N 



+ .9(1 - a) 
+ .9(1 



(27) 







N - 1 



)fori ^ 1 



Hu = Hn = - — fori ^ 1 



Hi 



The two lowest energy eigenvalues are 



/ 90 
N {N-l) 



fori l,j 7^ 1 



= i Uf + g)± J(f~(a-0)g) 2 + ^(a-0)fg\ , 



(28) 



and the higher eigenvalue E n = f + g is (N — 2)-fold degenerate Note that the minimum energy gap between 
the lowest two states is: 

g mm = \E+ - E_\ min ~ 0(1/ VN) (29) 
which goes to zero for large N, and explains heuristically the scaling of the running time. 



F. Constant-Time-Deutsch-Jozsa Algorithm 

Since the standard Deutsch-Jozsa algorithm requires only one call of the oracle operator, it is natural to ask whether 
one can do better than a vN improvement over the classical case. This question has been answered in the affirmative 
in [2£| where the authors were able to construct an algorithm that solved the problem in constant time, independent 
of N. In order to achieve this, they start with the same initial state, but choose the final state to have the form 

N/2-1 N/2-1 

Ilk) = -?== Y |2z> H 7 £= Y \2i + l) (30) 

This version of Deutsch's algorithm is optimized in the sense that the running time can be shown to be independent 
of N. The Hamiltonian is given by: 

1 N ~ X 
{j',fc}=o 
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2a N/2 ~ 1 2/3 N/2 ^ 

{Hi) mn — S mn — — S m (^2j)S n (2k) — -JT ^m(2j + l)^n(2fc+l) 

b",fe}=o {i:fc}=0 

In this case, a straightfoward calculation reveals the two lowest energy eigenvalues to be: 

E ± = \±^\-2s(\~ S ) (32) 

This yields a minimum energy gap that is independent of N which at a heuristic level accounts for the fact that the 
running time is independent of N . 



III. QUANTIFYING THE POSSIBLE RESOURCES: ENTANGLEMENT AND FIDELITY 

A. Entropy of entanglement 

There is no generally accepted measure for the entanglement of an arbitrary state. The exception is the case of a 
bipartite system in a pure state. 

For a bipartite system, composed of subsystems A and B, the Hilbert space is given by a tensor product Ha <S> Hb- 
The state of the composite system is a product state if it can be written as a product of pure states of subsystems, 

\^ab) = \*a) ® (33) 

If the state of the system cannot be written as a product state, it is entangled. An example of an entangled state for 
the simplest bipartite system, a system of two qubits, is the singlet Bell state 

l*>Bell - -^{\°U ® |1>B - |1>A ® \0)b) ■ (34) 

The von Neumann entropy as a measure of entanglement makes use of the fact that for an entangled pure state, 
the state of the composite system is known completely, but the states of the subsystems are not (they are not in pure 
states) . It measures the amount of information about one subsystem that can be obtained by making a measurement 
on the other subsystem. It is calculated from the reduced density matrix of either of the subsystems: 

SvN(Pred) = -Tr (predlog 2 Pred) (35) 

= -^// n l0g 2 /in, 

n 

where the fj, n are the eigenvalues of the reduced density matrix. If the composite state is a product state, it has zero 
entanglement and the reduced density matrix will have eigenvalues zero and one. From Eq. ll.'iot the von Neumann 
entropy will be zero. If the composite system is not in a product state, the von Neumann entropy will be non-zero. 
Its maximum value for a bipartite system is 1. 

If the system is in a pure state, but is multipartite, then it is not clear how it should be partitioned. In other 
words, the von Neumann entropy can depend on which qubits are traced out. It is unclear what physical significance 
to attach to an entanglement measure that depends on the specific partitioning of the system. For a multipartite 
system, or for mixed states, several other, inequivalent measures of entanglement have been suggested. Although there 
are general conditions that need to be satisfied by any viable measure (see, e.g., for a discussion), no generally 
accepted unique definition exists as yet. Several authors have considered many-particle systems, e.g. spin-systems, as 
bipartite by simply splitting them in the middle Another possibility is to consider one qubit in the "mean field" 
generated by the other qubits, as it is done, for example, in |2£l |. 

In this paper we study systems with n = 2 and n = 3 qubits. We give explicit expressions for the eigenvalues of the 
reduced density matrix for n = 2 qubits and n = 3 qubits in the Appendix B. 



B. Measuring entanglement geometrically 

In this section we suggest a different function to measure the entanglement of a state. This definition essentially 
measures the distance between the composite state under consideration and the product state that is closest, in a 
space parametrized by the coefficients of this product state. We note that this is a physically motivated approach, 
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different from the mathematical or axiomatic [27j, or the operational approach (as suggested in |29j). One advantage 
of the definition that we propose in this section is that its generalization for the case of more than two subsystems is 
straightforward: the problem of how to partition the system does not occur. For simplicity of presentation we discuss 
only real cofficients. 

We begin by discussing the Hilbert-Schmidt measure of entanglement Dhs, which is defined as the minimum 
(over the set S of product states) of the Hilbert-Schmidt distance between the density matrix p c of the state under 
consideration and the density matrix p p of a product state 30] : 

Dhs(Pc) = m.md H s{p P ,Pc), (36) 
where the Hilbert-Schmidt distance is defined as 



d H s{Pp,Pc) = \\p P - Pc|| ; \\p\\ = s/Trp^p. (37) 
To simplify the notation, we discuss below the case n = 2. A general product state of two qubits can be written: 

\ P ) = (ai|0) + a 2 |l))O(6i|0) + 6 2 |l)) (38) 
= a x 6i|00) + ai6 2 |01) + a 2 6i|10) + a 2 6 2 |ll) 
= (ai6i,ai6 2 ,a 2 6i,a 2 6 2 ) 

(Note that throughout this paper we write the components of vectors as rows of the form (v%, v 2 , V3, • • •), instead of 
as columns, to save space). The function to be minimized is then (the square root of) 

d 2 HS {ai, a 2 , 61, b 2 ) = 2 - 2(ai6ic + aib 2 ci + a 2 b 1 c 2 + a 2 b 2 c 3 ) 2 . (39) 

Using this definition as a starting point, we define a new function as follows. The difference between the state under 
consideration, |c), and a product state is given by: 

\d) = |c) - \p) = (c - aih, ci - a x b 2 ,c 2 - a 2 b^c z - a 2 b 2 ) (40) 

Our entanglement measure is the minimum of the distance (d\d) for all product states. The function to be minimized 
is: 

£(01,02,61,62) = (d\d) (41) 
= c 2 + c\ + c\ + c\ + [a\ + a 2 2 ){b\ + b\) 
— 2{a\b\CQ + aib 2 ci + a 2 b\c 2 + a 2 b 2 cz), 

where the state under consideration is taken to be normalized, which gives c§ + c 2 + c\ + c 2 . = 1. The minimum of 
this function with respect to the parameters of the product state is a measure of how much the state |c) differs from 
a product state, and thus a measure of entanglement. One advantage of our definition D in comparison with the 
Hilbert-Schmidt distance is that it is faster to calculate numerically, basically because of the absence of crossterms 
(compare Eqs. 1(39)1 and lUDl l. The advantage over the von Neumann entropy is that there is no ambiguity about 
how to partition a multi-partite system. 

The next step is to minimize the function 1141)) with respect to the parameters of the product state. We discuss this 
procedure below for the case n = 2. In Appendix C we consider n — 3. We will consider two different procedures. 
The first method uses normalized product states only. This approach is motivated by the fact that physical states 
are normalized. The normalization condition for the product state is 

(oifei) 2 + (ai6 2 ) 2 + (o 2 &!) 2 + (a 2 b 2 ) 2 = N a ■ N b = 1 (42) 

with N a = a\ +a 2 , Nb — b 2 + b 2 . We do not assume that the one-qubit states are individually normalized. For the case 
of n — 2 we have obtained analytic expressions for the values of the parameters {ai, a 2 , 61, 6 2 } that correspond to the 
product state closest to any arbitrary state |c). These expressions are given in Appendix D. We have also obtained 
numerical solutions for n — 3. These results are shown in section l)IV!l . 

We also study the minimum of l)lT)l using unnormalized product states. In this case it is easy to show that, in 
order find a non-trivial minimum of the distance function, one must satisfy a condition on the normalization constants 
N a ,Nb. We obtain this condition by looking at the four equations 

dD__dD__dD__dD_ {] 
dai da 2 db\ db 2 
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and rearranging terms. The result is that, in order to find a non-trivial solution, we must satisfy the condition: 

(A^ Q iV b ) 2 -^ a iV b + | CoC 3- ClC2 | 2 = (44) 

We note that this equation has the same form as the equation from which we obtain the eigenvalues of the reduced 
density matrix (see Appendix B.l): the equation det(p2 ed — = nas the form 

p 2 - p + |coc 3 - cic 2 | 2 = 0. (45) 

This result suggests that for n = 2 there is a connection between entanglement as defined through the von Neumann 
entropy, and entanglement defined as the minimum distance to an unnormalized product state. To understand this 
connection, it is useful to consider a geometric interpretation (see Appendix D for details). The space of all normalized 
4-dimensional states |c) is a 3-sphere. The space of normalized product states is a subset on the surface of this 3- 
sphere. Unnormalized product states lie along radial lines that intersect the 3-sphere. The closest product state to an 
arbitrary normalized state will lie on such a radial line in the interior of the 3-sphere ie. it will be an unnormalized 
product state. The only case where the closest product state will be normalized is if the original arbitrary state is 
itself a product state. This follows directly from Eq. 114411 . which implies that N a N(, < 1 unless \c0C3 — c\C2\ 2 = 0. 
Since the quantity N a N\, is the length of the 4-dimensional product state vector and it satisfies the same equation as 
the eigenvalues of the reduced density matrix, En. lEll provides a geometric interpretation of one of these eigenvalues. 
Therefore, one of the terms of the von Neumann entropy for an arbitrary normalized state is determined by the length 
of the closest non-normalized product state. The equations that are analogous to <44ll for n = 3 are found in Appendix 
C. 

We obtain numerical solutions for D m i„ as a function of time, using both normalized and unnormalized product 
states. In all cases, the results are almost identical. We show both graphs for the case of the n — 2 search algorithm 
only. 



C. Fidelity 

The fidelity measures the closeness of two states, in the sense defined below. The general definition of the fidelity 
of two states is written in terms of their density matrices p and a as [3ll ] 



F(p,a)=tr^^pa^p. (46) 
For pure states p = and a = the fidelity is just the overlap 

F(\i)Ao)) = Mj)\ (47) 

of the two states. In our case, we are interested in measuring the closeness of the instantaneous eigenstate \c(t)) and 
the final state |c(T)). We look at the function 

2V-1 



F(t) = 1 - F(\c(t)), |c(T)» = 1 - £ c*{t)ci{T) (48) 



i=0 



which approaches zero as t — > T. Introducing a slight abuse of terminology, we henceforth refer to F as the fidelity. 
Using this convention, the fidelity approaches zero as the desired state is reached (i.e. in the limit t — > T). 



D. Another distance measure 



There are several other natural definitions for the distance between the instantaneous ground state \c(t)} and the 
final state \c(T)}. We propose the following function: 



2V-1 

G(i) = £11^)1-1^)11 

i=0 



(49) 



10 



IV. CALCULATIONS AND RESULTS 

We work with a set of initial states that have zero initial entanglement. Our goal is to study the relationship 
between the running time of the algorithm and the amount of entanglement that is produced during the evolution of 
the algorithm. As described in section i|IIBJl . the initial state determines the Hamiltonian H(s) for each algorithm as 
a function of the dimensionless time variable s — t/T. The eigenvectors and eigenvalues can be obtained numerically 
at each point in time. From these expressions we can extract the coefficients }co, ci, ■ ■ -cjv_i} that correspond to 
the components of the eigenvectors in the computational basis (see Eqns l|A4l) and ijAfiH . Using these coefficients we 
calculate the von Neumann entropy from l-'>5ll , (IB7II , l|B10ll ; the geometric entanglement as the minimum of l|39|) ; the 
fidelity from I|48J1 : and the geometric separation from l|49|i . These quantities are obtained numerically as functions of 
the time s throughout the evolution of the system from the initial state \ipo) at s = to the final state at s = 1. 

We also present some results in terms of the optimized time t. As explained in section ijllCJl . the optimized time is 
obtained by applying the adiabatic condition 11 fill to each infinitesimal time step dt throughout the evolution, instead 
of globally to the entire time interval T. The function s(t) = t/T is replaced with a non-linear function that is 
calculated numerically. 

We note that the search algorithm and Deutsch's algorithm with {a = 1, (3 = 0) are identical since we can choose 
(without loss of generality) the marked state of the search algorithm to be \m) = |0). As a consequence, we only 
have to investigate the case (a = 0, (3 = 1) for Deutsch's algorithm. We also note that for the search algorithm and 
the constant time Deutsch algorithm, the system evolves to a final state which has zero entanglement (since the final 
state is a product state). On the other hand, Deutsch's algorithm with (3 = 1 evolves to a final state with non-zero 
entanglement. 

The initial states that we use, for n = 2 and n = 3, are given below: 

n = 2 (50) 

|red) = Normalized{l, 3/2, 1, 3/2} 

| yellow) = Normalized{l, 1,4/3,4/3} 

| green) = Normalized} 1 , 1,1,1} 

|blue) = Normalized}!, 1, 2/3, 2/3} 

| cyan) = Normalized} 2, 2,1,1} 

|magenta) = Normalized}3, 1,3, 1} 

n = 3 (51) 

|red) = Normalized} 1 , 1, 3, 3, 1, 1, 3, 3} 

|yellow) = Normalized}3, 3, 3, 3, 1,1,1, 1} 

|green) = Normalized} 1 , 1, 1, 1, 1, 1, 1, 1} 

|blue) = Normalized} 1, 1, 2, 2, 1, 1, 2, 2} 

|cyan) = Normalized}2, 2, 2, 2, 3, 3, 3, 3} 

(magenta) = Normalized} 19, 19, 1, 1, 19, 19, 1, 1} 

The general features of our results are as follows: 

[A] Consistency checks on the calculation: 

(1) The optimized running time is less than the unoptimized running time. 

(2) The running time is larger where the initial fidelity is larger. 

(3) The same behaviour is seen for the geometric separation. 

[B] Search Algorithm: 

(1) Using the unoptimized time, the region where the magnitude of the rate of change of the fidelity is greatest 
corresponds to the region of greatest entropy. 

(2) Using the optimized time, the speed of the algorithm is adjusted so that it moves faster at either end of the time 
interval and more slowly in the middle, where the gap is smallest. This effectively flattens out the curve that gives 
the rate of change of the fidelity. 

(3) Using optimized time, states with a larger maximum entanglement tend towards longer running times. We 
remind the reader that the interpretation of this result is unclear, because of the inherent structure of the adiabatic 
calculation: different initial states have different initial fidelity, and thus it is not clear how to separate the effects of 
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fidelity and entanglement on the running time. In other words, it seems likely that the longer running time is actually 
caused by the larger initial fidelity, and not by the larger maximum entanglement. This conclusion is supported by 
point (B.l) which indicates that increased entanglement causes the fidelity to move more quickly towards its final value. 

[C] Deutsch Algorithm: 

(1) Using the optimized time variable we find that states that generate entropy more quickly have a shorter running 
time. 



[D] Constant Time Deutsch Algorithm: 

(1) Using the unoptimized time, the region where the magnitude of the rate of change of the fidelity is greatest 
corresponds to the region of greatest entropy. 

(2) Using optimized time, states with a larger maximum entanglement tend towards longer running times. As 
discussed for point (B-3) above, it seems likely that the longer running time is actually caused by the larger initial 
fidelity, and not by the larger maximum entanglement. This conclusion is supported by point (D.l). 

We give below the numerical results and graphs that illustrate these points. 



A. Search Algorithm for two qubits 



The running times are (point (A-l)): 



T 



unoptimizcd 



optimized 



{Tred = 598, Tyollow = 503, Tgreen = 346, Tbluo = 234, T cyan = 194, T r 



magenta 



{^rcd — 234, Ty e ii ow — 213, Tgreen — 173, Tb me — 137, Tcyan — 122, Tj 



magenta 




PIG. 1: Largest initial fidelity corresponds to the largest running time (point (A-2)). 



1 . 5 



0.5 




50 100 150 200 



FIG. 2: Largest initial geometric separation corresponds to the largest running time (point (A-3)) 
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FIG. 3: The region of greatest magnitude of the rate of change of the fidelity corresponds to the region of greatest entanglement 
(point (B-l)). This graph is for the red initial state. 




t-red 



PIG. 4: The rate of change of the fidelity as a function of the optimized time is flattened (point (B-2)). This graph is for the 
red initial state. 



SvN (t) 




t 



50 100 150 200 



FIG. 5: Entropy versus optimized time (point (B-3)). 



B. Search Algorithm for three qubits 

The running times are (point (A-l)): 




100 200 300 400 500 600 



FIG. 8: Largest initial fidelity corresponds to the largest running time (point (A-2)). 



G(t) 




100 200 300 400 500 600 



FIG. 9: Largest initial geometric separation corresponds to the largest running time (point (A-3)). 
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FIG. 10: The region of greatest magnitude of the rate of change of the fidelity corresponds to the region of greatest entanglement 
(point (B-l)). This graph is for the red initial state. 




PIG. 11: The rate of change of the fidelity as a function of the optimized time is flattened (point (B-2)). This graph is for the 
red initial state. 



SvN (t) 




t 



100 200 300 400 500 600 



FIG. 12: von Neumann entropy versus optimized time (point (B-3)). 



C. Deutsch Algorithm for two qubits 



running times are (point (A-l)): 




— 55, -Zred — 

= 49, T red = 



FIG. 14: Largest initial fidelity corresponds to the largest running time (point (A-2)). 



G(t) 



1 . 4 
1.2 
1 




FIG. 15: Largest initial geometric separation corresponds to the largest running time (point (A-3)). 




FIG. 17: States that generate normalized geometric entanglement more quickly have a shorter running time (point (C-l)). 



D. Deutsch Algorithm for three qubits 

running times are (point (A-l)): 



T. 



unoptimizcd 



optimized 



{^magenta — 217, Xyellow — 104, 2~red — 55, T\>\ uc — 43, Tg^een — 40, Xcyan — 
{-^magenta 131, '7y C [} ow 81, T rc d 49, ?bi uc 39, Tg rccn 37, T C y an 
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FIG. 18: Largest initial fidelity corresponds to the largest running time (point (A-2)). 



G(t) 
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PIG. 19: Largest initial geometric separation corresponds to the largest running time (point (A-3)). 




FIG. 21: States that generate geometric entanglement more quickly have a shorter running time (point (C-l)). 



E. Constant-Time-Deutsch Algorithm for two qubits 

The running times are (point (A-l)): 




FIG. 22: Largest initial fidelity corresponds to the largest running time (point (A-2)). 



G(t) 



1 . 4 

1.2 
1 




FIG. 23: Largest initial geometric separation corresponds to the largest running time (point (A-3)). 
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S(s) , F' (5) 



-0 . 1 
-0.2 
-0.3 
-0 . 4 
-0.5 



FIG. 24: The region of greatest magnitude of the rate of change of the fidelity corresponds to the region of greatest entanglement 
(point (D-l)). This graph is for the cyan initial state. 



SvN (t) 




FIG. 25: von Neumann entropy versus optimized time (point (D-2)). 



V. DISCUSSION AND OUTLOOK 

The purpose of this paper is to study the role of entanglement in quantum computing. We have used the adiabatic 
formulation of quantum computing, so that we can explicitly track the time evolution of the quantities we are interested 
in. We have studied two different algorithms, adiabatic quantum search and the adiabatic Deutsch-Jozsa algorithm. 
We have also considered a variation of the Deutsch-Jozsa algorithm j23| which we call the constant-time-Deutsch- 
Jozsa algorithm. We have introduced a new function to describe the physical property of entanglement and compared 
with the results for the entropy of entanglement. We also introduced a different function characterizing the distance 
between states and compared the results with the traditional fidelity. We have calculated the evolution of these 
quantities as a function of time, and also as a function of an optimized time variable, which corresponds to applying 
the adiabatic approximation in a maximally efficient way (T^ |. 

For the search algorithm, we find that the region where the magnitude of the rate of change of the fidelity is greatest 
corresponds to the region of greatest entropy production. This result indicates that the state approaches the final 
state most quickly when the entanglement is greatest, which supports the idea that entanglement is a resource for 
speed-up in quantum computation. We note that there also seems to be a correlation between greater maximum 
entanglement and longer running times. Although this result appears to contradict the conclusion above, we believe 
that it is an artifact of the structure of the adiabatic computation. We expect that both the initial fidelity, and the 
amount of entanglement that is produced during the evolution, will affect the running time of the algorithm. However, 
since the Hamiltonian, and therefore the algorithm, depends explicitly on the choice of initial state, we are unable to 
study the roles of these factors independently. In other words, we believe that the longer running times produced in 
the cases we have studied are the result of larger initial fidelity, and not driven by larger maximum entanglement. 

The results for the constant-time-Deutsch-Jozsa algorithm are similar to the results for the search algorithm. 
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The Deutsch-Jozsa algorithm is different from the previous two cases in that the entanglement of the final state is 
not zero. In this case, using the optimized time variable, we find that the states that generate entropy more quickly 
have a shorter running time. This result also seems to support the idea that entanglement is a resource for speed-up 
in quantum computation. 

In this paper we have worked only with n = 2 and n — 3. The next step is to extend our calculations to systems 
with larger n. Work on this is in progress. 



APPENDIX A: NOTATION AND COMPUTATIONAL BASIS 



The basis states for one qubit in the S z representation are 

|0) = 

|1) = 



(Al) 



where we use the symbol = to denote a representation. In the computational basis, an n-qubit state is written as 

\x) n = \x xi...x n ) (A2) 

where Xi = or 1 V i, and the string they form is the binary representation of x. We explicitly write the index n 
here to denote the number of qubits. The dimension of the Hilbert space for a system of n qubits is N = 2". 



The computational basis for n — 2 qubits is: 



|0>2 


= 100) 


= |0>i <s 


5|0)i 


= (1,0,0,0) 


|1>2 


= |01) 


= |0)i«5 


5|l)i 


= (0,1,0,0) 


|2>2 


= |10) 


= |l>i€ 


5|0)i 


= (0,0,1,0) 


|3>2 


= 1") 


= |l>i€ 




= (0,0,0,1) 



A general state of two qubits is written as 



|*>2 = ^2ci\i) 2 = (Cq, Ci, c 2 , c 3 ) 



i=0 



The computational basis for n = 3 qubits is: 



0) 3 = 


|000) = 


|0)i€ 


) |0)i <i 


5|0)i 


= (1,0,0,0,0,0,0,0) 


1)3 = 


|001) = 


|0)i€ 


) |o)i a 


5|l)i 


= (0,1,0,0,0,0,0,0) 


2)3 = 


|010) = 


|0)i€ 


) |i>i <§ 


H0)i 


= (0,0,1,0,0,0,0,0) 


3)3 = 


|011) = 


|0)i€ 


) |l)i<5 


Mi)i 


= (0,0,0,1,0,0,0,0) 


4) 3 = 


|100) = 


|l)i€ 


) |0)i§ 


Mo)i 


= (0,0,0,0,1,0,0,0) 


5)3 = 


|101) = 


|l)i€ 


) |0)i§ 


Mi)i 


= (0,0,0,0,0,1,0,0) 


6)3 = 


|110) = 


|l)i€ 


) |iM 


Mo)i 


= (0,0,0,0,0,0,1,0) 


7)3 = 


1111) = 


|l)i€ 


) |l>i <5 


Mi)i 


= (0,0,0,0,0,0,0,1) 



A general state of three qubits is written as 

7 

|*)3 = y^Cj|j) 3 = (c , Ci, C 2 , C 3 , C 4 , C 5 , C 6 , C 7 ) 



(A3) 



(A4) 



(A5) 



(A6) 



»=o 



APPENDIX B: DENSITY MATRICES AND REDUCED DENSITY MATRICES 



The density matrix for an ensemble of pure states {pi,\ipi)} is 

P = ^2pi\ipi)(ipi 



(Bl) 
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meaning that the system is in state with probability pi. Recall that Trp = J^iPi = 1 an d Trp 2 > 1. Since p is 
positive, it must have a spectral decomposition 

p=X>ii>oi> ( B2 ) 

where the vectors |j) are orthogonal, and Xj are real, non-negative eigenvalues of p (cf., e.g., 31j) . 

For the special case of a pure state (or ensemble), one of the Pi is unity and all others are zero. The density matrix 
reduces to 

pP ure =|^i)<^|. (B3) 

In this case only, Tr(p purc ) 2 = 1. The eigenvalues of the pure state density matrix are and 1. 

The reduced density matrix is used for the description of subsystems of a composite quantum system. Assume we 
have a composite system with the subsystems A and B. It is described by the density matrix p AB . The reduced 
density operator for system A is defined as 

p A = Tr B {p AB ) 1 (B4) 
where Ttb is the partial trace over systems B (cf.,e.g.,[3l))- 



1. Calculations for n — 2 qubits 

We construct the density matrix for a pure state of two qubits (i.e. a composite system) using 

EH and £b3J|. The 

density matrix is 



Pi 



= |*>22<*|= $>SWI 



1,1 



( C Cq CqC* C cl CqcI 

CiCq CiC* CiC* 2 C1C3 

C 2 Cq C2C\ C2C2 C2C3 

V C 3 C5 C 3 C* C3C2 C3C3 



(B5) 



We obtain the reduced density matrix for one of the subsystems by tracing out the degrees of freedom corresponding 
to the other qubit: 



red _j_ I C C 

92 ~ 1 c 2 c* 



CiC-l C C 2 
C 3 C* C2C2 



ClC 3 
C3C3 



(B6) 



This reduced density matrix has the eigenvalues 

1 



/i± = -(1 ± v/l-4|c C3-CiC 2 | 2 ) 



(B7) 



with p + + p- = 1. Note that for the case of a two qubit system, the reduced density matrices obtained by tracing 
out either of the two subsystems are equal. 



2. Calculations for n — 3 qubits 



For a pure composite state of three qubits we use <HU) and lO to get: 



P3 = |*)3 3<*|=$^CiC;|*) 3 3<7l 





c c\ 


C0C2 


C0C3 


C0C4 


c ct, 




C C? 


\ 




ClC* 




C1C3 


C1C4 


C1C5 




C1C7 




c 2 c* Q 




C2C2 


C2C3 


C2C| 


C2C 5 


C2C6 


c 2 c 7 




C3C5 


C3C1 


C3C2 


C3C3 


C3CI 


* 

C3C5 


C3C 6 


C3C7 




C A Cq 


C 4 C* 


C4C2 


C4C3 


C4C4 


c 4 c; 


c 4 c; 


C4Cy 




C5C0 


C 5 C* 


C5C2 


C5C3 


C5C4 


C5C5 




C5C7 






C 6 C* 


C 6 C2 


CqC* 3 


C6C4 


C6C5 




C6C7 


/ 


V C 7C(5 


C 7 C* 


C7C2 


C7C3 


C7C4 


C 7C 5 


C7Cq 


C7C7 



(B8) 
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For the three qubit system we obtain the reduced density matrix by tracing out the degrees of freedom corresponding 
to the second and third qubits. We obtain: 



The eigenvalues are: 



red ^ / C C + C X C X + C 2 C 2 + C3C3 C C 4 + C x C b + C 2 C 6 + C 3 C 7 , 
1 ' ' C4Cq + C 5 C* + C 6 C 2 + C 7 C* 3 CiC% + C5C5 + C 6 Cg + C7C7 ' ' " 



1 



/i± = -\l±y/l- AXXj ; M+ +^_ = l (BIO) 

XX = (Fx • V?)(V 2 ■ V 2 *) - (Vt • V 2 *)(V 2 • t?) 
Vi = (c , ci, c 2 , c 3 ) 

V 2 = (C4,C 5 ,C 6 ,C 7 ) 

(Note that this formula also works for the two qubit case with the notation V\ — (c , c x ) and V 2 = (c 2 , c 3 ).) 
If we trace out the first and second qubit, the reduced density matrix becomes 

-red ± f CqCq + C 2 C 2 + C±C% + C§C% C Q c\ + C 2 Cg + C4C5 + C G C* 7 \ 

3 \ ciCq + c 3 c 2 + c 5 c* 4 + c 7 c% cxc\ + C3C3 + C5C5 + c 7 c£ ) 

In this case, the eigenvalues are given by Eq. ijBlOl! as well, but with 

Vi = {c ,c 2 ,Ci,c§) (B12) 
V 2 = (d,c 3 ,c 5 ,c 7 ). 

If we trace out the first and third qubit, the reduced density matrix becomes 

-red ^ / C Cq + C X C\ + C±c\ + C 5 C^ C C 2 + CiC* 3 + C 4 Cg + C 5 C 7 \ . . 

\ C 2 C*q + C 3 C* + C 6 cl + C 7 C% C 2 C 2 + C3C3 + C 6 Cg + c 7 c 7 J \ ' 



and the eigenvalues are given by Eq. IIBIOJ) with 

V 1 = (co,ci,c 4 ,c 5 ) (B14) 

^2 = (C2,C 3 ,C 6 ,C 7 ). 

APPENDIX C: CONSISTENCY CONDITION FOR AN UNNORMALIZED THREE QUBIT PRODUCT 

STATE 

The general expression for a three qubit product state is given by (compare (l38l) l: 

\ P ) = (ai|0) + aa|l))®(6i|0) + 62|l))®(ei|0) + ea|l)) (CI) 
= {ai&id, ai&ie 2 , a\b 2 et, atb 2 e 2 , a 2 biei, a 2 bte 2 , a 2 b 2 ei, a 2 b 2 e 2 } 

We do not assume the product state is normalized. The geometric distance between this state and an arbitrary state 
(co, ci, C2, c 3 , C4, C5, C6, C7) is (compare (@H): 

D = (aibiei - c ) 2 + (ai&ie 2 - ci) 2 + (ai& 2 ei - c 2 ) 2 + (ai& 2 e 2 - c 3 ) 2 (C2) 
+ (a2&iei - c 4 ) 2 + (a 2 bie 2 — c 5 ) 2 + {a 2 b 2 ei — c 6 ) 2 + (a 2 b 2 e 2 — c 7 ) 2 

We can minimize this expression with respect to the parameters (at, a 2 , b\, b 2 , ei, e 2 ). Doing so results in a set of 6 
equations, for which one can find three consistency equations by demanding that non-trivial solutions exist. Defining 

N a =a\+a 2 2 , N b = b\ + b 2 2l N e = e\ + e 2 , (C3) 

these equations are 

= N A e N 2 a N% (C4) 
- N 2 N a N b [(cl + c% + cl + cl)e\ + ((c x c Q + c 7 c 6 + c 3 c 2 + c 5 c 4 )eie 2 ) + (c\ + c 2 + c 2 + c 2 )e 2 ] 
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+ ((c4C 2 - c c 6 )e 2 + (c 3 c 4 + c 5 c 2 - c 6 ci - c 7 c )eie 2 + (c 5 c 3 - cic 7 )ef) 

= A^iN^e 

- N 2 N b N e [(c 2 + c 2 2 + c 2 3 + cg)a? + ((c 4 c + c 6 c 2 + c 5 ci + c 3 c 7 )aia 2 ) + (c 2 - + c 2 + c 2 + c 2 6 )a 2 2 ] 
+ ((-C2C1 + c 3 Co)a 2 + (-cic 6 + c c 7 + c 3 c 4 - c 2 c 5 )aia 2 + (c 4 c 7 - cec^a 2 ) 2 

= <iV 2 7V 2 

- N 2 N a N e [{cl + c 2 + c\ + cl)b\ + ((c 7 c 5 + c 2 c + c 3 ci + c 4 c 6 )6i& 2 ) + b\ + (c 2 + c 2 + c 2 + c 2 )6 2 ] 
+ ((-C4C1 + c 5 c )6 2 + (-c 6 ci + c c 7 - c 4 c 3 + c 2 c 5 )6i6 2 + (c 2 c 7 - c 3 c 6 )bl) 2 

The relationship of these consistency conditions to the different ways of tracing out states when defining the von 
Neumann entropy is currently under investigation. 



APPENDIX D: ANALYTIC SOLUTIONS FOR MINIMUM DISTANCE TO A NORMALIZED TWO 

QUBIT PRODUCT STATE 

For n — 2 the general product state is given in ffiSjl. If we use the change of variables: 

ai = ^f~N~ a cos9 a 2 = y/N^sin9 ( D1 ) 
b\ = \f~N~b cos 4> hi = \fNbS'm(f> 

where: 

N a = a\ + a% and N b = b\ + bl. 

the product state becomes: 

\J N a N b (cos cos <p, cos 9 sin <p, sin 9 cos 0, sin 9 sin </>) . (D2) 

The direction of the state 4-vector is fully determined by 9 and (j>. The normalizations A^ a and Ni, of the individual 
states determine the length of the state 4-vector. The normalization condition for the product state is N a Nb = 1. 
The distance between a general normalized state: 

{c , ci, c 2 , c 3 } where c 3 = 1 - c 2 a - c 2 - c 2 (D3) 

and a general product state is: 

(d\d) = (cq - aifoi) 2 + (ci - ai& 2 ) 2 + (c 2 - a 2 foi) 2 + (c 3 - a 2 & 2 ) 2 . (D4) 
Using the change of variables above leads to: 

(d\d) = (c - ^N a N b cos 9 cos 2 + (c x - ^N a N b cos 9 sin A 2 (D5) 
c 2 — y/ N a N b sin cos (/)] + (c 3 — \J N a N b sin 9 sin 



In order to find the minimum distance to a product state we want to find values for 9 and (j> such that the derivative 
of this distance with respect to both variables (9 and (j>) vanishes. Differentiating this equation with respect to 9 and 
setting the derivative equal to zero leads to the equation: 

tan^ C2 + C3tan t (D6) 
Co + c 4 tan (p 

Differentiating this equation with respect to 4> and setting the derivative equal to zero leads to the equation: 

+ x ci+c 3 tan0 /-n^ 

tan^= -. (D7) 

Co + c 2 tan 9 

Neither of these equations depends on the normalization N a N b of the general product state. Solution of these equations 
therefore leads to a line of product states with arbitrary length and the intersection of this line with the unit 4-sphere 
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is the physical product state. Note also that when the state under consideration is a product state, it will satisfy 
c c 3 = C1C2, and in this case llPfil) and IID7JI reduce to: 

tan#=— ; tan0= — (D8) 
Co c 

which is consistent with cq = a\b\, c\ = a\bi and C2 = ct2&i, and the definitions in l|Dl|l , 
Substituting (IPfSf) into l|D7|) and solving for tan0 leads to the solutions: 



,1+4- (cl+cl) ± V / 4(coc 1 +c 2C3 ) 2 + (c2 + c|-(c2 + c |)) 2 



tan</> = y —— . (D9) 

2 (c ci + c 2 c 3 ) 



For an arbitrary normalized state determined by Co, ci and C2, the closest product state can be determined by ijDfiJ) 
and (IDOI . The distance to this product state can then be determined from ijDfiJl after choosing its normalization, 
N a N b . 



APPENDIX E: EQUATIONS FOR MINIMUM DISTANCE TO A NORMALIZED THREE QUBIT 

PRODUCT STATE 

For n — 3 the general product state is given in ijClJl . If we use the change of variables: 

ai = \fW a cos9 a 2 = \fW a sm9 (El) 

bi = \f~N~b cos 4> bi — \/~N b sm.(p 

ei = \J N e cos ip 62 = \J ~N~ e sin ip 

where: 



N a = a\ + a\ N b = b\ + b\ N e = e j + e\. 



the product state becomes: 



\/ N a N b N e (cos 9 cos (p cos ip, cos 9 cos <f> sin ip, cos 9 sin (p cos ip, cos 9 sin (p sin ip, (E2) 
sin 9 cos 4> cos ip, sin 9 cos (p sin -0, sin 9 sin <p cos ip sin 9 sin <p sin ip) . 

The direction of the state 8-vector is fully determined by 9, <p and ip. The normalizations N a , Nb and N e of the 
individual states determine the length of the state 8-vector. The normalization condition for the product state is 

N a N b N e = 1. 

The distance between a general normalized state and a product state is given in l|C2Jl . Using the change of variables 
above leads to: 

(d\d) = (c - ^N a N b N e cos 9 cos <p cos ip) + (c ± - ^N a N b N e cos 9 cos (p sin ip) (E3) 

+ ( C2 — \J N a N b N e cos 9 sin <p cos + (ca - y/N a N b N e cos 9 sin cp sin ip ) . 

+ (c 4 - y/N a N b N e sin 9 cos cos sin cos sin ^ • 

+ ^C6 — \/ N a N b N e sin 6* sin cos + — \J N a N b N e sin 9 sin sin ipj . 

In order to find the minimum distance to a product state we want to find values for 9, <p and ip such that the derivative 
of this distance with respect to these variables (6, (p and ip) vanishes. Differentiating with respect to each variable 
and setting the derivative equal to zero leads to the equations: 

C4 + cq tan 6 + C5 tan ip + C7 tan <b tan ip 
tan 9 = ■ ■ ■ . (E4) 

co + C2 tan (p + ci tan ip + C3 tan <fi tan ip 

C2 + C6 tan 6* + C3 tan ip + c-j tan 9 tan ip 

tan0= ■ -. (E5) 

Co + C4 tan 9 + c\ tan ip + C5 tan ft tan t/j 
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ci + C5 tan 9 + C3 tan 6 + cj tan 6 tan 6 . 

tan?/; = ■ -. (E6) 

Co + C4 tan + C2 tan <p + c§ tan tan 

None of these equations depends on the normalization N a NbN e of the general product state. Solution of these 
equations therefore leads to a line of product states with arbitrary length and the intersection of this line with the 
unit 8-sphere is the physical product state. Solutions of these equations can be substituted into l|E3Jl to find the 
distance to the nearest product state. 

This set of equations is easily generalizable to the case for arbitrary n. It is not clear that analytic solutions of 
these equations would be possible for arbitrary n, but they can always be solved numerically. They can be used to 
provide a measure of the entanglement which is symmmetric among the qubits. The advantage of this approach is 
that it provides an unambigous measure of entanglement. 
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